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Abstract 

A Higgs-like particle with a mass of about 126 GeV has been discovered at the LHC. Within 
the experimental uncertainties, the measured properties of this new state are compatible with those 
of the Higgs boson in the Standard Model (SM). While not statistically significant at present, the 
results show some interesting patterns of deviations from the SM predictions, in particular a higher 
rate in the 77 decay mode observed by ATLAS and CMS, and a somewhat smaller rate in the 
t + t - mode. The LHC discovery is also compatible with the predictions of the Higgs sector of the 
Minimal Supersymmetric Standard Model (MSSM), interpreting the new state as either the light 
or the heavy CP-even MSSM Higgs boson. Within the framework of the MSSM with seven free 
parameters (pMSSM-7), we fit the various rates of cross section times branching ratio as measured 
by the LHC and Tevatron experiments under the hypotheses of either the light or the heavy CP-even 
Higgs boson being the new state around 126 GeV, with and without the inclusion of further low- 
energy observables. We find an overall good quality of the fits, with the best fit points exhibiting 
an enhancement of the 77 rate, as well as a small suppression of the bb and r + r _ channels with 
respect to their SM expectations, depending on the details of the fit. For the fits including the 
whole data set the light CP-even Higgs interpretation in the MSSM results in a higher relative fit 
probability than the SM fit. On the other hand, we find that the present data also permit the 
more exotic interpretation in terms of the heavy CP-even MSSM Higgs, which could give rise to 
experimental signatures of additional Higgs states in the near future. 
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1 Introduction 



One of the most important goals of the Large Hadron Collider (LHC) is to explore the origin of 
electroweak symmetry breaking (EWSB). The spectacular discovery of a Higgs-like particle with a 
mass around Mh — 126 GeV, which has been announced by ATLAS [I] and CMS |2j, marks a 
milestone of an effort that has been ongoing for almost half a century and opens up a new era of 
particle physics. Both ATLAS and CMS reported a clear excess in the two photon channel, as well as 
in the ZZ^ channel. The discovery is further corroborated, though not with high significance, by the 
VFVF(*) channel and by the final Tevatron results [3]. The combined sensitivity in each of the LHC 
experiments reaches more than 5 a. The observed rate in the 77 channel turns out to be considerably 
above the expectation for a Standard Model (SM) Higgs both for ATLAS [4] and CMS [5], whereas 
the bb and the t + t - channels appear to be somewhat low in the LHC measurements [6,7J. While 
those possible deviations from the SM prediction are not statistically significant at present, so that 
within the uncertainties the results are compatible with the SM, if confirmed in the future the observed 
patterns could be a first indication of a non-SM nature of the new state. 

Among the most studied candidates for EWSB in the literature are the Higgs mechanism within 
the SM and within the Minimal Supersymmetric Standard Model (MSSM). Contrary to the SM, two 
Higgs doublets are required in the MSSM, resulting in five physical Higgs boson degrees of freedom. 
Without explicit CP- violation in the soft supersymmetry-breaking terms these are the light and heavy 
CP-even Higgs bosons, h and i7, the CP-odd Higgs boson, A 1 and the charged Higgs boson, H ± . The 
Higgs sector of the MSSM can be specified at lowest order in terms of the Z boson mass, Mz, the 
CP-odd Higgs boson mass, Ma, and tan/3 = i^/^i, the ratio of the two Higgs vacuum expectation 
values. The masses of the CP-even neutral Higgs bosons and the charged Higgs boson can be predicted, 
including higher-order corrections, in terms of the other MSSM parameters, see, e.g., [8>9| for reviews. 

It was shown that in particular the interpretation of the new state as the light CP-even Higgs boson 
of the MSSM is a viable possibility (called the "light Higgs case" in the following). The implications 
and phenomenology of this scenario has been studied in a series of papers [TUl - fH] . On the other hand, 
it was also pointed out that the heavy CP-even Higgs boson can have a mass around 126 GeV [TOlfTT] 
(called the "heavy Higgs case") while maintaining (within the uncertainties) a SM-like behaviour. All 
five MSSM Higgs bosons in this scenario would be rather light, and it would in particular imply the 
existence of another light Higgs with a mass below 126 GeV and suppressed couplings to W and Z 
bosons. For a recent discussion of the phenomenology of such a scenario, see also [15] . 

The question arises whether the MSSM (or another model beyond the SM) can give a prediction 
of the production cross sections and decay widths of the observed Higgs-like state that yields a better 
description of the data than the one provided by the SM. While at the current level of accuracy no 
clear deviation from the SM can be claimed, the situation could change once the whole dataset of 
~ 20 fb _1 to be collected in 2012 is incorporated in the analyses. The main aim of this paper is to 
investigate in how much the MSSM can improve the theoretical description of the current experimental 
data, and potentially which parts of the parameter space of the MSSM are favoured by the current 
experimental data from the various Higgs search channels. 

Because of the large number of free parameters, the MSSM Higgs search results at LEP [16], 
the Tevatron [17] and the LHC [18J have been interpreted in certain benchmark scenarios [T9l - [2T] (of 
which the m^ 13 ^ scenario has been most widely used), where the MSSM parameters entering via higher- 
order contributions are set to specific fixed benchmark values. However, in order to investigate the 
potentially favoured regions in the MSSM parameter space a scan over the relevant SUSY parameters 
has to be performed. A complete scan over the in principle more than hundred free parameters of the 
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MSSM parameter space is neither technically feasible nor does the available experimental information 
provide sufficient sensitivity to simultaneously constrain a large number of parameters. One therefore 
needs to focus on a certain subset of parameters. The most ambitious scans performed up to now were 
done [22] for the phenomenological MSSM with 19 free parameters (pMSSM-19, see [23J for details). 
However, on the one hand it is difficult to sample such a multi-dimensional parameter space sufficiently 
well, on the other hand it is well known that several of the parameters of the pMSSM-19 hardly affect 
Higgs phenomenology. We therefore focus in this paper on a smaller set of parameters, namely the 
phenomenological MSSM with the seven free parameters that we regard as most relevant for the 
phenomenology of Higgs and flavour physics (pMSSM-7, see below for details on these parameters). 
This seven-dimensional parameter space, which as we will demonstrate captures most of the allowed 
Higgs phenomenology of the MSSM, can be sampled quite well with O(10 7 ) scan points. We will 
comment in our analysis on the potential impact of varying further MSSM parameters. 

In our analysis we perform fits in the MSSM both for the interpretation of the LHC signal in 
terms of the light and the heavy CP-even Higgs of the MSSM and we compare the fit results with 
the SM case. We take into account all available individual search channels at ATLAS and CMS 
at 7 and 8 TeV center-of-mass energy that have been published by the end of July 2012, including 
also the two combined Higgs mass values. We furthermore include the final results of the Tevatron, 
corresponding to three additional channels. Besides the direct Higgs search channels we take into 
account the Higgs exclusion bounds, limits on the SUSY masses, as well as the most relevant set of 
low -energy observables, BR(6 — y ^t)? BR(i? s — y fjL~^ fi ), BR(5 U — y tz/ t ), — 2)^ and the mass of the 
W boson, Mw- 

The paper is organised as follows: Section 2 gives a summary of the most relevant supersymmetric 
sectors and parameters. In section 3 we briefly review the calculations and codes used for the Higgs 
sector predictions, as well as for the low-energy observables. Details on the parameter scan are 
given. In section 4 the results of the scan, the best-fit points and the preferred parameter regions are 
presented. We briefly discuss the origin of the most significant deviations of the pMSSM-7 predictions 
with respect to the SM results. The conclusions can be found in section 5. 

2 Theoretical framework for pMSSM-7 

In the following we briefly describe the relevant sectors of the MSSM and define the seven basic 
parameters for our scan. As a first and general simplification we restrict ourselves to the MSSM with 
real parameters. The restriction to seven free/ independent parameters is based on the intention to 
sample the full parameter space with (9(10 7 ) scan points, while ensuring that the most important 
effects in the MSSM Higgs sector and flavour physics are covered. 

The tree-level values for the predictions of the MSSM Higgs sector quantities are determined by 
tan/3, the CP-odd Higgs-boson mass Ma, and the Z boson mass Mz- These predictions include the 
other Higgs boson masses, their couplings to other MSSM particles, their production cross sections, 
as well as their decay properties. Consequently, we choose as free parameters the two tree-level 
parameters, 

(i) : M A , 

(ii) : tan/3. 

Beyond tree-level, the main correction to the Higgs boson masses stems from the t/t sector, and 
for large values of tan/3 also from the b/b sector. In order to fix our notations, we list the conventions 
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for the input parameters from the scalar top and scalar bottom sector of the MSSM. The notation 
can be taken over to scalar taus and scalar tau neutrinos via the substitutions t —± z> T , b —> f, taking 
into account that the MSSM contains only the scalar superpartners of the left-handed neutrinos. 
Furthermore, the same notation holds for the first and the second generation of scalar fermions. The 
mass matrices in the basis of the current eigenstates £l,£r and are given by 



KA 2 I M h +™? + cos2/3(±-§4)Mf mtX t 
Mi ~ \ m t X t M^ + m t 2 + |cos2/34Ml I' 



U 2 I M ! L + m l + cos 2/3(-i + |4) M l m "X b \ 

M b ~ y m b X b Ml R + ml-\ cos 2/?4Mf J' {2) 

where M^, Mj R , and M^ r are the soft SUSY-breaking mass parameters in the scalar top and 



bottom sector, m t and 7715 are the respective quark masses, s w — y 1 — M^/M^ with Mw denoting 
the mass of the W boson, and 

m t X t — m t (A t — fi cot f3) , — (A^ — fi tan /3) . (3) 

Here A t (A^) denotes the trilinear Higgs-stop (Higgs-sbottom) coupling. The higgsino mass parameter 
/i, which also appears in Eq. ([3D, is taken as free scan parameter 

(hi) : /x. 

SU(2) gauge invariance leads to the relation M$ L = . Using a universal parameter also for the 
left /right-handed squark soft masses of the third generation (and similarly for the sleptons), we choose 
as free parameters 

(iv) : M~ q , := M~ tL (= M-J = M 1r = M^, 

(v) : M l3 :=M fL (=M, L ) = M fR1 

(vi) : Af := A t = A b = A T . 

For the soft scalar masses of the first two generations, which are much less relevant for Higgs physics 
(but can play a relevant role for the low-energy observables) , we choose fixed parameter values as 

Mq L = Mq R (q = c, 5, u, d) = 1000 GeV, (4) 
M Jl = M Ir (I = fi : Vp, e, u e ) = 300 GeV. (5) 

The choice for the first and second generation squarks places their masses roughly at the level currently 
probed at the LHC. Somewhat larger values would have a minor impact on our analysis. The values 
for the first and second generation slepton mass parameters were chosen to provide rough agreement 
with the anomalous magnetic moment of the muon (see below). 

The trilinear Higgs coupling parameter for the first two generations we set to 

A c ,s,u,c,iJ,,e — Af. (6) 
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Parameter Minimum Maximum 



M A [GeV] 90 1000 

tan/3 1 60 

fji [GeV] 200 3000 

Mq 3 [GeV] 200 1500 

M~ h [GeV] 200 1500 

A f [GeV] -3Mq 3 3Mq 3 

M 2 [GeV] 200 500 



Table 1: Ranges used for the free parameters in the pMSSM-7 scan. 



The Higgs sector observables furthermore depend on the SU(2) gaugino mass parameter, 



(vii) : M 2 



which we take as the final free parameter in our analysis. The other electroweak gaugino mass 
parameter, Mi, is fixed via the GUT relation 



The gluino mass parameter, which enters the Higgs mass predictions only from two-loop order on, is 
fixed to a value close to the limits from recent searches at the LHC, 



An adjustment of the gluino mass parameter to even higher values is expected to have a negligible 
impact on our analysis. 

3 Model predictions and constraints 
3.1 Parameter space scanning 

The pMSSM-7 parameter space is sampled by performing random scans (using uniform distributions) 
over the seven input parameters in the ranges given in Tab. [TJ The two cases, where either h or 
H corresponds to the signal at Mh ~ 126 GeV, are treated in two separate scans, and the results 
are discussed in parallel below. Each scan starts with O(10 7 ) randomly chosen points with a flat 
distribution over the parameter ranges. Dedicated, smaller, sampling is then performed to map the 
interesting regions of parameter space0 In practice, the full parameter ranges from Tab. [1] are taken 
only for the light Higgs case, while for the heavy Higgs case we limit Ma < 200 GeV and tan /3 < 30 
(still using the full ranges for the other parameters), which improves the sampling efficiency in the 
relevant mass region for Mh- In addition to the free parameters listed in Tab.[U the remaining 
parameters are fixed according to section [2l 

Since we are mainly interested in the Higgs sector, we do not exploit the full possibilities in the 
low-energy MSSM to vary the soft-breaking parameters of the first two generations or the gluino mass 
(we will comment below on a possible impact of this choice on the MSSM predictions for (g — 2) M ). 

lr The reader should keep in mind here (and in the following) that the point density has no statistical meaning. 




(7) 



M 3 = m~ g = 1000 GeV. 



(8) 
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Consequently, it is not relevant to apply LHC exclusion bounds from supersymmetry searches in 
our analysis, since these can be avoided by adjusting the additional parameters to sufficiently high 
values with only small effects on the Higgs sector. We do, however, apply the results from Higgs 
boson searches, see the next subsection. We also apply the model-independent limits on sfermion and 
chargino masses, typically at the level of ~ 100 GeV from direct searches at LEP (as summarized 
in the PDG review [24]). Furthermore, we require that the lightest supersymmetric particle (LSP) 
is the lightest neutralino. The top quark pole mass is sampled from a Gaussian distribution with 
m t = 173.2 ± 0.9 GeV, using a cutoff at ±2 a. Effects of other parametric uncertainties from SM 
quantities are estimated to be small, and are therefore neglected. 

For the evaluation of the sparticle and Higgs masses we use the code FeynHiggs (version 2.9.4) [25]- 
[27] . The residual Higgs mass uncertainty from this calculation (e.g. from missing higher orders) is 
estimated to be around 2 GeV, depending on the considered region of parameter space [26J. We are 
interested in parameter points that give a Higgs mass prediction, for either or Mh , close to the 
observed LHC signal. We therefore constrain the analysis in a first step to points with or Mh in the 
region 121 — 129 GeV. To avoid configurations in parameter space that give an unstable perturbative 
behavior in the Higgs mass calculation, we use a criterion based on the Z-matrix (as defined in [27J) 
and exclude points for which \ \Z^\ — \Zj^ || /\Zj^\ > 0.1 (see [TT] for a similar treatment). Here k — 1 
(2) is set for a SM-like light (heavy) Higgs, and the superscripts refer to the same quantity evaluated 
with 1-loop (1L) or 2-loop (2L) precision. 

To obtain an indication of what the currently favoured regions of the MSSM parameter space are, 
we use a simple statistical treatment of the data where the different observables are taken into account 
by calculating, for every parameter point in the scan, a global x 2 function 



X 2 = £ + ^ + + ^ ( 9 ) 

=1 1 1=1 1 M H i=l 1 



(/iz-A,) 2 , y- v ~ AO 2 , (M h , H - M H f ' ^° (gj - o % 



Quantities with a hat denote experimental measurements, and unhatted quantities the corresponding 
model predictions for Higgs signal strength modifiers, and low-energy observables (LEO), 0{. The 
different observables entering Eq. ([9j) are described in more detail in the following sections. The 
combined uncertainties G{ contain the known theory and experimental uncertainties. Correlations are 
neglected throughout, since they are for most cases not publicly available. The total number of degrees 
of freedom, z/, is counted in the naive way as v — n G b s — %>ara 5 where n G b s = ^lhc + ^Tev + 1 + ^leo 
(for LHC, Tevatron, the Higgs boson mass, and low-energy observables); n para is the number of model 
parameters. In the SM we have n para = 1 (the Higgs mass), and for both MSSM analyses n para = 7. 



3.2 Direct Higgs searches 

We use two kinds of experimental data from direct Higgs searches in our analysis: exclusion limits at 
95% confidence level (CL), and the signal rates measured in different channels for the observed LHC 
signal around Mh ~ 126 GeV. The exclusion limits from LEP, the Tevatron, and the LHC (including 
data published before ICHEP 2012) are taken into account using HiggsBounds version 3.8.0 [28] . Since 
no x 2 information is available from these searches, but only the exclusion limits at a fixed level of 
confidence, we apply these limits as hard cuts on the parameter space. It should be noted, however, 
that HiggsBounds only tests for each parameter point the model predictions against the single channel 
with the highest expected sensitivity for an exclusion, in order to ensure a consistent interpretation 
of the exclusion limit as a 95% CL. 
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The measured signal strength modifiers, for the observed Higgs-like state around Mh ~ 
126 GeV are taken into account in our fit directly in the x 2 evaluation (see Eq. ([9j)). The data 
for all included channels is given in Tab. EJ with the corresponding experimental signal strengths and 
their (asymmetric) 1 a error bands. These rates provide the main dataset to which we fit the MSSM 
Higgs sector. In total we include 37 observables, where 34 are from the LHC experiments and 3 pro- 
vide supplementary information from the Tevatron. The best fit signal strength modifiers of ATLAS 
and CMS are given for different Higgs masses, corresponding to the values measured by the individual 
experiments, i.e. we interpret the experimental discoveries as being compatible, and due to a single 
new state. The Tevatron data, which does not admit a mass measurement from the observed excess 
on its own, is evaluated for Mh = 125 GeV. All values listed in Tab. [2] are extracted directly from the 
quoted experimental references, with one exception: ATLAS has not provided a measurement for the 
signal strength modifier of H —> ZZ^ separately for the 7 and 8 TeV data, but only for the combina- 
tion (the 7 TeV values are available from a previous analysis). To compare to our 8 TeV predictions, 
these values are therefore calculated from the 7 TeV and 7 + 8 TeV data under the assumption of 
independent Gaussian measurements, following the procedure outlined in [29]. This should lead to an 
uncertainty on the estimated 8 TeV rate of the same order as the overall uncertainty from neglecting 
the (unknown) correlations. 

The MSSM predictions for the signal strength modifiers are evaluated according to 

Efe Uik(?k(pp -+h,H) x BR(/i, H -> i) 

where (Jk(pp —> h, H) denotes the contribution to the Higgs production cross section from partonic 
subprocess fc, evaluated at the predicted Higgs mass. The production modes considered are gluon- 
gluon fusion (gg), vector boson fusion (VBF), associated vector boson production (VH), and associated 
tih(H) production. The experimental efficiencies cj^ have only been published by ATLAS and CMS 
for the 77 analysis; for CMS in the case of the subcategories, and for ATLAS also the inclusive result. 
We make use of these numbers when they are available. For all other channels we have to use the 
"naive" efficiencies deducible from the analysis description (e.g. for a VBF-type analysis tagging two 
forward jets, we set uj — 1 for the VBF cross section, whereas all other modes have cj = 0). In 
channels where the mass resolution is not good enough to separate contributions from different Higgs 
bosons, we approximate the contributions from H and the CP-odd Higgs A by adding their signal 
rates incoherently. We do not add the rates of the CP-even Higgs bosons, whose joint contributions 
to the signal could also include interference effects. Our analysis is therefore limited to the case with 
a single CP-even Higgs boson close to the observed signal, and we leave a more detailed treatment 
of the case when ~ Mh ~ 126 GeV to a dedicated analysis. Since the CP-odd Higgs does not 
have tree- level couplings to vector bosons (and hence also a reduced coupling to photons), it gives 
a negligible contribution to the channels with vector bosons in the Higgs production and/or decay. 
Effectively, the CP-odd Higgs therefore only plays a role for the inclusive (0/1 jet) t + t~ channels. 
In these channels it can easily dominate over the H contribution for large values of tan f3. In the 
light Higgs case, we find that the masses of h and A differ by Ma — M^ > 50 GeV in the favoured 
region (see below). Thus we do not take any contributions to the h rates from the CP-odd Higgs into 
account. 

The cross section predictions entering Eq. ([TO]) are calculated, both in the MSSM and the SM, 
using FeynHiggs (version 2.9.4) [25H2T] . For the SM cross sections the results of the LHC Higgs 
cross section working group are implemented [34] (where the gg production cross sections are taken 
from [35j). The corresponding MSSM production cross-sections are obtained in the effective-coupling 
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Channel 


[TeV] Alow 






Reference 


ATLAS data at M H -- 


= 126.5 GeV 




bb 


7 -1.646 


0.510 


2.680 


m 


TT 


7 -1.550 


0.464 


2.011 


n 


WW 


7 -0.164 


0.438 


1.103 


L30J 


WW 


8 1.308 


1.920 


2.536 


m\ 


77 (inclusive) 


7 1.397 


2.155 


2.903 


m 


77 (inclusive) 


8 1.054 


1.685 


2.326 


m 


ZZ 


7 0.405 


1.080 


2.177 


[32j 


zz 


8 0.400 


1.049 


1.708 


133] 




CMS data at M H = 


125.0 GeV 




bb (VH) 


7 -0.606 


0.588 


1.824 


m 


bb (VH) 


8 -0.441 


0.424 


1.535 


17J 


bb (till) 


7 -2.624 


-0.771 


1.288 


[7j 


tt (0/1 jet) 


7 -0.400 


1.000 


2.441 


[7j 


tt (0/1 jet) 


8 0.588 


2.153 


3.635 


[7j 


tt (VBF) 


7 -2.912 


-1.718 


-0.359 


[7J 


tt (VBF) 


8 -3.035 


-1.759 


-0.400 


[7J 


TT (VH) 


7 -2.418 


0.671 


4.788 


[7J 


77 (Dijet loose) 


8 -2.660 


-0.626 


1.409 


[5j 


77 (Dijet Tight) 


8 -0.267 


1.289 


2.868 


15J 


77 (Untagged 3) 


8 2.007 


3.754 


5.549 


[5J 


77 (Untagged 2) 


8 -0.195 


0.930 


2.080 


15J 


77 (Untagged 1) 


8 0.475 


1.504 


2.533 


15J 


77 (Untagged 0) 


8 0.212 


1.456 


2.701 


15J 


77 (Dijet) 


7 2.174 


4.209 


6.243 


[5j 


77 (Untagged 3) 


7 -0.099 


1.528 


3.132 


15J 


77 (Untagged 2) 


7 -0.434 


0.715 


1.887 


15J 


77 (Untagged 1) 


7 -0.291 


0.643 


1.600 


[5J 


77 (Untagged 0) 


7 1.337 


3.132 


4.974 


15J 


WW (0/1 jet) 


7 -0.029 


0.588 


1.206 


[7J 


WW (0/1 jet) 


8 0.176 


0.835 


1.494 


[7J 


WW (VBF) 


7 -3.900 


-1.306 


0.918 


[7J 


WW (VBF) 


8 -0.523 


1.371 


3.347 


[7J 


WW (VH) 


7 -5.753 


-2.829 


0.341 


[7J 


ZZ 


7 0.176 


0.671 


1.371 


[7J 


ZZ 


8 0.259 


0.794 


1.494 





Tevatron data at Mjj 


= 125.0 GeV 




bb 


1.96 1.290 


1.970 


2.710 


131 


77 


1.96 1.080 


3.620 


6.580 


L3J 


WW 


1.96 0.000 


0.320 


1.450 


131 



Table 2: Experimentally measured values for the Higgs signal strength modifiers, and their corre- 
sponding uncertainties (lower/upper edges of la error bars), in the various channels. 
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approximation [35]. The gg production cross section follows the description in [37], where results 
of [38] were used. The decay width evaluation includes a full one-loop correction for the decay to 
fermions [391S0] ; see [37] for more details on the other channels. 

In addition to the signal strength modifiers, we include a x 2 contribution from the measured Higgs 
mass Mh (which in the MSSM can correspond to either Mh or M/j). Unlike in the SM, Mh and 
Mn are not free parameters in the MSSM but a prediction of the theory. Averaging the ATLAS and 
CMS measurements, we obtain Mjj — 125.7 GeV. To fully cover the difference between the individual 
experimental measurements we assign a conservative uncertainty of = 1 GeV. With more data 
we expect this experimental uncertainty on the Higgs mass to be significantly reduced. In addition, 
we add linearly to this number a theoretical uncertainty cr^ eo = 2 GeV [26J, which accounts for the 
uncertainty in the MSSM Higgs mass calculation. The total uncertainty entering Eq. ([9j) is therefore 

a M H = 3 GeV ' 

3.3 Low-energy observables 

In addition to the measurements related to the LHC Higgs signal, we include several low-energy ob- 
servables (LEO) in the fit. These are listed in Tab. El which summarizes the experimental valueqj (for 
the case of BR(£> S the upper limit) and the corresponding SM theory predictions (evaluated 

for Mj^ M = 125.7 GeV and m t = 173.2 GeV). The flavour physics observables are evaluated (both in 
the SM and the MSSM) using Superlso (version 3.2) [32], which in particular contains the results 
for BR(£> — »• X s j) based on the NNLO calculation of [43J. Our fit includes also the anomalous mag- 
netic moment of the muon, = \{g — 2)^, which shows a deviation of more than 3a between the 
experimental measurement and the SM prediction!! We use Superlso to calculate the MSSM contri- 
bution da^ to the anomalous magnetic moment of the muon (we have cross-checked those results with 
FeynHiggs and found good agreement) including the dominant two-loop contributions [45] , see [46] 
for a review. 

As a final observable we also include the MSSM prediction of the W boson mass into the fit. Here 
the SM prediction shows a ~ 1.5 a deviation from the latest experimental value, M^ 9 — 80.385 ± 
0.015 GeV [47]. Our MSSM evaluation of M w is done using FeynHiggs, where the full SM result [38] 
is supplemented with the leading corrections from the t/b sector [9 , 49 j . A comparison with the best 
available MSSM evaluation [50] shows that corrections larger than 10 MeV can be missed if some 
uncoloured SUSY particles are light. Consequently, we assign a theory uncertainty of 15 MeV to our 
Mw evaluation and conservatively combine it with the experimental uncertainty linearly. Thus in 
total we take an uncertainty of ±30 MeV into account. 



2 We note that the Belle Collaboration has recently reported a new measurement of BR(B n — >• ru T ) that is in better 
agreement with the SM (and thus naturally with models with two Higgs doublets) [41 . While we do not take this new 
result into account in our overall fit results, we do comment briefly on its possible effects. 

3 The most recent evaluation of a M in the SM [44], taking also r data into account, finds an even larger deviation of 
more than 4 a. 
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Observable 


Experimental value 




SM value 


BR(B -> X s7 ) 


(3.43 ± 0.21 ± 0.07) x 10" 4 


m 


(3.08 ± 0.22) x 10" 4 


BR(B S — »> a* + A»~) 


< 4.2 x 10" 9 


M 


(3.55 ± 0.38) x 10" 9 


BR(£ U -> r^ T ) 


(1.66 ± 0.33) x 10" 4 


M 


(1.01 ±0.29) x 10" 4 




(30.2 ± 9.0) x 10" 10 


[SUED 




il% 


(80.385 ± 0.015) GeV 


[47j 


(80.363 ± 0.004) GeV 



Table 3: The experimental values and (SM) theory predictions for low-energy observables (LEO) used 
to constrain the MSSM parameter space. 



Case 


LHC only 

x 2 l v xl P 


LHC+Tevatron 

x 2 l v xl P 


LHC+LEO 

x 2 l v xl P 


LHC+Tevatron+LEO 

x 2 l v xl P 


SM 
h 
H 


27.6/34 0.81 0.77 
23.3/28 0.83 0.72 
26.0/28 0.93 0.57 


31.0/37 0.84 0.74 
26.8/31 0.86 0.68 
33.1/31 1.07 0.37 


41.6/39 1.07 0.36 
26.7/33 0.81 0.77 
35.5/33 1.08 0.35 


45.3/42 1.08 0.34 
30.4/36 0.84 0.73 
42.4/36 1.18 0.21 



Table 4: Global y 2 results with v degrees of freedom from the fits of the SM and the MSSM with 
either h or H as the LHC signal, the reduced xl = X 2 / U i and the corresponding p- values. The number 
of degrees of freedom are evaluated naively as v — n G b s — %>aram- 



4 Results 

In Tab. S] we present the results of our fits in terms of total x 2 values (with v degrees of freedom), the 
reduced xl = X 2 l v -> an d the corresponding p- values. Since v is derived via the naive counting, the 
absolute numbers of the p-values should not be overinterpreted; the relative numbers, however, give 
a good impression of the relative goodness of the fits. For each MSSM intrepretation (the cases of 
either h or H as the 125.7 GeV signal) we present four different fits: one taking the complete dataset 
(LHC+Tevatron+LEO) into account, one where the low-energy observables (LEO) are left out, one 
where the Tevatron data are left out, and finally the fit where only LHC observables are considered. 
When the fit is performed using only the high-energy collider data, both with and without the Tevatron 
results, the obtained x 2 values of the best fit points are quite similar between the SM and the two 
MSSM interpretations, where the fit in the heavy Higgs case becomes slightly worse after the inclusion 
of the Tevatron data. When low energy observables are included, the SM and the heavy Higgs case 
fits become somewhat worse. In the latter case this can be understood from the potentially larger 
contributions of light Higgs bosons to £>-physics observables. For the SM fit the reason lies in the fact 
that the SM prediction for (g — 2)^ differs by more than 3 a from the experimental value. Still we find 
that the SM provides a good fit to the full dataset, with psm — 0.34. On the other hand, concerning 
the MSSM it should be kept in mind that we did not fit the second generation slepton masses, which 
could potentially further improve the a M fit. For the complete fit, the corresponding p- values in the 
MSSM cases are ph = 0.73 (pn = 0.21) for the h (H) interpretations, respectively, which are both 
acceptable p-values. Overall, the data shows no clear preference for the MSSM over the SM at this 
point. While the MSSM fits, in particular in the light Higgs case, yield lower x 2 values than the SM, 
this comes at the expense of additional parameters, so that the difference in the p-values is rather 
moderate. It is interesting to note that the fit to the heavy Higgs case, while not being as good as 
in the light Higgs case, still provides an acceptable description of all the data, i.e. the MSSM in the 
non-decoupling region still provides a viable solution. This interesting case will be analyzed in more 
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Figure 1: Fit results for the signal strength modifiers, in the case that the light CP-even Higgs 
is interpreted as the new boson around ~ 125.7 GeV ("light Higgs case"). The experimental data is 
shown as black dots (with error bars). The other symbols show best fit points, corresponding to the 
full fit (LHC+Tevatron+LEO) (blue solid squares), without the Tevatron data (blue open squares), 
and without LEO (red diamonds). 
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Figure 2: Fit results for the signal strength modifiers, in the case that the heavy CP-even Higgs 
is interpreted as the new boson around ~ 125.7 GeV ("heavy Higgs case"). The experimental data is 
shown as black dots (with error bars). The other symbols show best fit points, corresponding to the 
full fit (LHC+Tevatron+LEO) (blue solid squares), without the Tevatron data (blue open squares), 
and without LEO (red diamonds). 
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detail below. 

Starting with the best fit for the h case, we show in Fig. []] the different best fit points using 
all available data (LHC, Tevatron, LEO) (blue solid squares), leaving out LEO (red diamonds) or 
leaving out the Tevatron data (blue open squares). The comparison of these three different types of 
results allows to trace the origin of trends in the fitted parameters. The experimental data on the signal 
strength modifiers in the different channels (as indicated in the figure) is shown as black dots, with the 
error bars corresponding to ±1 a uncertainties on jl. The values for the best fit point of the complete 
fit (LHC, Tevatron, LEO) are also presented in tabular form in Tab. El From here we can determine 
some characteristics of the best fit point, such as a significantly enhanced rate in the 77 final state 
and nearly SM rates for the other channels. Leaving out the Tevatron data a (small) suppression of 
the fermionic final states can be observed. The fitted rates demonstrate that the pMSSM-7 is able to 
accomodate the main trends in the LHC/Tevatron data. Comparing the best fit points with/without 
LEO, we find a qualitatively very similar behaviour. 

In Tab. [5] we give the details on the results for the low-energy observables. In the light Higgs case, 
the only relevant contribution to the total x 2 comes from BR(B U — »• ru r ). Using an experimental value 
close to the new Belle result [41] would substantially reduce this x 2 contribution and lead to an even 
better fit. The best-fit value of BR(£> S —> /i + /i~) lies below the SM prediction. This feature is indeed 
found for most of our favoured region. We have checked that this trend is present already without 
taking the x 2 contribution of BR(£> S — > /i + itself into account, see also the discussion in [14] . 

The best fit points for the heavy Higgs case are presented in Fig. [2] (numerical values in Tab. Ej). 
As the figure shows, the same best fit point (albeit with different total x 2 ) is obtained for the different 
cases with/without LEO. Leaving out the Tevatron data, however, has a larger qualitative impact on 
the results, and rates close to zero are allowed in the bb channel, which will be discussed below in more 
detail. Since the alternative hypothesis of a heavy Higgs explanation for the LHC signal is fitted to 
the same data as for the h case, and we have already seen that the overall x 2 is similar to (although 
slightly higher than) the light Higgs case, it is perhaps not so surprising to find that comparable rates 
are obtained for the best fit in the heavy Higgs case. This illustrates again that this interpretation 
is also a possible scenario to explain the LHC Higgs signal in the MSSM. The more in-depth results 
below are therefore presented in parallel for the two separate cases with h or H corresponding to the 
signal discovered at the LHC. 

In Tab. [5] we also give the results for the low-energy observables in the heavy Higgs case. One 
can see that the relatively small value of the Higgs mass scale in this case leads to non-negligible x 2 
contributions from BR(£> —> X.7) and BK(B U —> rv T \ where the latter would substantially improve 
for a value close to the new Belle result. Also the SUSY contribution to a M turns out to be relatively 
small, giving a sizable contribution to the total x 2 (which is however affected by our choice to keep 
the slepton mass parameters fixed). Concerning BR(£> S /i + /x~) it should be noted that, as in the 
light Higgs case, the preferred value is below the SM result, which again holds for most of the favoured 
region. 

We now turn from the global fit properties and the best fit points to a more detailed analysis of 
the scan results. Fig. [3] shows distributions of Ax\ = Xh ~ Xh m in (light Higgs case) for the different 
signal rates. The colour coding is as follows: all points analyzed in the scan (which pass theoretical 
consistency checks and have one CP-even Higgs boson in the interval 121 GeV < Mh < 129 GeV) are 
shown in gray. The blue points in addition fulfill constraints at 95% CL from direct Higgs searches 
applied by HiggsBounds 3.8.0. The three rows show the results for the full fit (upper row), excluding 
LEO (middle row), and excluding the Tevatron data (lower row). The signal rates are calculated as the 
inclusive Higgs production cross section (evaluated at = 8 TeV) times the decay rate, normalized 
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Table 5: Best fit results (for the complete fit) with corresponding x 2 contributions and pulls for each 
observable. 
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Figure 3: Distributions of Ax\ versus the different signal rates (defined in the text) for the light 
Higgs case. The colours show all points in the scan (gray), and points that pass the direct Higgs 
search constraints from HiggsBounds (v. 3.8.0) (blue). The three rows show the distributions (with 
the corresponding minimal x 2 subtracted) for the complete fit (upper row), excluding LEO (middle 
row), and excluding Tevatron data (lower row). 



to the SM predictions 



K x - 



h,H) x BR(/i, H -> X) 



°i(pp ^h,H)x BR (/i, H -> X) 



(11) 



The only final state for which we consider a different observable than the fully inclusive Higgs produc- 
tion is 66, where the sum is only taken over the cross sections for (h,H)Z and (h,H)W ± associated 
production. As described above, for the inclusive t + t~ channels we consider the contribution of both 
H and A when these are close in mass. To make it clear when this is the case, we denote the joint (inclu- 

sive) rate as i? rr . We also define a common rate for vector boson final states Rvv '-— Rww — Rzz- 
To keep things simple, we do not include the experimental efficiencies for the 77 channel in Eq. (fTTI) . 
since the efficiencies are different for the two experiments. These are however used for the different 
predictions entering the fit (as described in Eq. (fT0|) ). Fig. Sallows to investigate the best fit rates 
in some more detail (subject to the approximations already discussed). Uncertainty intervals can be 
extracted from the range with < 1 (corresponding to 68% confidence intervals in the Gaussian 
case). The results for the complete fit are 
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(12) 
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Figure 4: Ax# versus the different signal rates (defined in the text) for the heavy Higgs case. Colour 
coding the same as in Fig. El The three rows show the distributions (with the corresponding minimal 
X 2 subtracted) for the complete fit (upper row), excluding LEO (middle row), and excluding Tevatron 
data (lower row). 

For R^ r we observe a distribution which is very flat near the minimum. This indicates a low sensitivity 
in the fit to constraints from t + t~ final states, and it permits substantially reduced t + t~ rates at 
a very low additional x 2 contribution. The light Higgs case could therefore easily explain strongly 
reduced t + t~ rates, although this is not visible in the best-fit point. 

Results for the heavy Higgs case are shown in Fig. [U The resulting Axjj distributions for individual 
Rx are similar to those for A%^, except for i? rr , where the additional contribution from the A boson 
strongly enhances this quantity over the corresponding result in the light Higgs case. Extracting the 
results for the minimal x 2 in the same way as for the light Higgs case, we obtain for the complete fit 

H H — 1 9^+0-30 r>H _ -i in +0.18 r>H _ n oo+0.05 U H I A — 9 p;/1+ - 31 (A<X\ 
J^VV — I'^-omi H ll - - LiU -0.06> K bb - U -^>-0.12> K tt -^O4_ 017 . (16) 

An interesting behaviour can be observed in the results leaving out the Tevatron data, as shown in 
the lower row of Fig. [U In that case the best fit points have a very small value for R^ h . More sampling 
would be needed to cover this particular case in detail. Since we focus in the following on results of 
the combined fit, we leave this for a dedicated study. 

More information about the phenomenology of the pMSSM-7 Higgs sector can be found from the 
correlations between the different rates. This is shown in Fig. [5] for the light Higgs case. Compared 
to the one-dimensional \ 2 distributions of Fig. El this figure introduces two new colours that are used 
in the following to show regions close to the minimum x 2 • We highlight points for which H < 2.3 
(red) and ^Xh h < ^.99 (yellow). In the Gaussian limit these correspond to 68% (95%) confidence 
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Figure 5: Correlations between signal rates for the light Higgs case. The colour coding follows that of 
Fig. El with the addition of the favoured regions with Ax 2 h < 2.3 (red) and Ax 2 h < 5.99 (yellow). The 
best fit point is indicated by a black star. 




R?y R b H b (VH) R b H b (VH) 



Figure 6: Correlations between signal rates in the heavy Higgs case. Colours similar to Fig. [5j but 
here representing Ax# < 2.3 (red) and < 5.99 (yellow). The black star indicates the best fit 

point for the heavy Higgs case. 



regions with two degrees of freedom. We shall refer to these points simply as the favoured region/points, 
or sometimes most favoured region/points when Ax\ H < 2.3 is discussed. Here (and in all figures 
from here on) we refer to the x 2 of the complete fit, including LHC, Tevatron and LEO. The best fit 
point is indicated in the figures by a black star. 

The left plot of Fig. [5] shows the strong, positive, correlation between Ry V and R^. In most of the 
viable parameter space we find R^ > Ry V . The favoured region contains points with fully correlated 
rates in the interval 0.9 < R^ lVV < 1.6, but also solutions with lower degree of correlation, where a 
77 enhancement (up to R^ ~ 1.8) is accompanied by a much smaller (or no) enhancement of Ry V . 
In the second plot of Fig. [5] we compare the results of i?^ 7 and R% b (we remind the reader that the 
latter rate is calculated using the VH production mode only). We find an anticorrelation between 
these two rates. This can be understood from the fact that the h,H —> bb decay gives the largest 
contribution to the total width for a Higgs boson in this mass range, both in the SM and (typically) 
also in the MSSM. A reduction of the h,H —> bb partial width is therefore effectively a reduction of 
the total decay width, which leads to a simultaneous enhancement of the branching ratios into the 
subdominant final states. This has recently been pointed out [TTJ[13] as an important mechanism to 
enhance the 77 rate in the MSSM. We shall see below how these effects on the Higgs decay widths 
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Figure 7: Effective couplings to vector bosons and gluon pairs of the Higgs boson corresponding to 
the LHC signal in the light Higgs case (left) and heavy Higgs case (right). The colour coding follows 
that of Fig. 

affect the parameters in our global fit. The third (right) plot in Fig. [5] shows the weak correlation 
of R^ T to R% h , where in principle any value of R^ T < 1 is found in the favoured region for R% b < 1. 
Consequently, it is possible to find a strong reduction of the t + t~ mode while maintaining a SM-like 
bb rate. 

Turning to Fig. [6j we show the rate correlations for the heavy Higgs case. Similar trends as in 
the light Higgs case are visible in the heavy Higgs data, with the notable difference in the t + t - rate, 
mainly due to the inclusion of the contribution from the CV-odd Higgs A. The favoured regions are 
found at values for RrJ A between 2 and 4, while R^ h remains below 1. 

Fig. [7| shows the effective couplings to vector bosons and gluon^l of the MSSM Higgs assigned 
to the LHC signal for both the light Higgs case (left plot) and heavy Higgs case (right plot). In 
both cases the favoured regions have effective Higgs couplings to vector bosons very close to the SM 
value (gfhHjyy — 1)? an d the light (heavy) CP-even Higgs boson behaves SM-like in this production 
(or decay) mode (which reflects the general fact that the VVh or VVH coupling cannot be larger 
than the corresponding SM coupling). For the effective coupling to gluons, g^ hH ^ gg , the somewhat 
larger favoured range 0.7 < dfh H)gg ~ ^'^ ^ s obtained. Consequently, an overall large enhancement 
of the Higgs production cross sections is disfavoured by the fit. The observed (positive) deviations 
from the SM in the VV and 77 channels for the best fit points can therefore be attributed mainly to 
changes in the decay branching ratios of the corresponding Higgs states, rather than modifications of 
the production rates. 

We now briefly discuss what mechanisms can alter the branching ratios in the manner observed, 
and what the consequences are for the favoured regions of MSSM parameter space. In Fig. [8] we show 
the scan results in the plane of the Higgs sector tree-level parameters (Ma, tan/?), where the results 
for the light (heavy) Higgs case are shown in the left (right) plot. In the light Higgs case one can 
note the region at low Ma, high tan/3, which is excluded by direct MSSM Higgs searches (mainly 
H/A — > t + t~). The excluded region appears smaller in this plane than the corresponding results 

4 The effective coupling to gluons is denned via the ratio of the MSSM and SM decay widths to gluons. The relative 
vector boson coupling is obtained from the inclusion of higher-order corrections into the CP-even Higgs mixing angle, a. 
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Figure 8: Higgs sector tree-level parameters (Ma, tan/3) in the light Higgs case (left) and in the heavy 
Higgs case (right). 

published by the experiments [18], since their results are shown only for one particular benchmark 
scenario (the so-called m™ ax scenario [20]). In an inclusive scan of the pMSSM-7 parameter space, 
points are found where higher order corrections to the bottom Yukawa coupling lead to suppressed 
production rates for the heavy MSSM Higgs bosons, and a larger fraction of the parameter space in 
the (Ma, tan /3) plane therefore opens up (see the analyses in [2HGS]). Sizable branching ratios of 
H/A to SUSY particles also reduce the sensitivity of the searches in the t + t~ final state. We see 
that the regions of very high tan/3 > 40, and also low tan/3 < 8, are disfavoured by the fit. At high 
tan /3 this results from a poor fit to (g — 2) M and flavour observables, whereas for low tan f3 the fit to 
the LHC Higgs observables becomes worse. For low tan (3 it also becomes increasingly difficult to fit 
the relatively high Higgs mass value (125.7 GeV), although viable solutions down to tan/3 ~ 4 can 
be found [10]. Low values of Ma are disfavoured by the fit results in the light Higgs case, with the 
preferred region starting at Ma > 170 GeV (and the most favoured region at Ma > 230 GeV). Taking 
the rate information into account therefore suggests somewhat higher mass scales for the MSSM Higgs 
sector than what is required by the M/j ~ 125.7 GeV Higgs mass measurement alone [10]. For the 
light Higgs case the lower limits on Ma in the favoured regions of the fit exclude the possibility of 
M H ± < m t , where the charged Higgs can be produced in the decay of the top quark. On the other 
hand, the region favoured by the fit does not show any upper limit for Ma, which demonstrates that 
the decoupling limit (corresponding to Ma ^> M^, where the MSSM Higgs sector reproduces the 
predictions for a SM Higgs) remains a possible scenario. This is to be expected, given the high quality 
of the SM fit to the LHC data. 

For the heavy Higgs case, as shown in the right plot of Fig. [HI the situation is very different. Low 
values for Ma are preferred, and the favoured region in (Ma, tan /3) is much smaller than for the light 
Higgs case: 110 GeV < M A < 140 GeV and 7 < tan /3 < 13. Even though the H can be very SM-like 
in this scenario, this situation is very different from the decoupling limit in the light Higgs case since it 
implies that all five MSSM Higgs bosons are light. In contrast to the light Higgs case, in this scenario 
values of the charged Higgs boson mass only below the top mass (M H ± < m t ) are found, which may 
offer good prospects for the searches for charged Higgs bosons in top quark decays. We therefore show 
in Fig. [9] the results for BR(t — )> bH + ) as a function of M H ±. The current upper limit on this decay 
mode is (9(1%) [56J, which is close to the maximal value favoured by the fit. With more integrated 
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M H+ (GeV) 

Figure 9: Branching ratio of the top quark into a charged Higgs boson and a bottom quark in the 
heavy Higgs case. 




M h (GeV) 

Figure 10: Effective coupling squared g\ zz of the lightest MSSM Higgs boson to a Z boson pair, as 
a function of the lightest Higgs mass in the heavy Higgs case (Mh ~ 125.7 GeV). 

luminosity, charged Higgs searches will therefore have an interesting sensitivity to probe the heavy 
Higgs scenario at the LHC in the near future. 

While in the heavy Higgs scenario the low preferred values for Ma typically lead to a situation 
where A, and are rather close in mass, the lightest Higgs boson, h, can have a significantly 
lower mass, as illustrated in Fig. [101 As we see from this figure, points with M^ < 90 GeV have a 
very small effective coupling to vector bosons, g\ zz <C 1, which explains why such light Higgs bosons 
are compatible with the Higgs search limits from LEP. The bulk of the favoured region is found for 
60 GeV < Mh < 90 GeV, with an effective coupling squared to vector bosons at the sub-percent level. 
Another feature which is clearly visible in the HiggsBounds allowed points (blue) is the degradation 
of the limit around M^ ~ 98 GeV, which was caused by a slight excess of events observed at LEP 
in that mass region. While a scenario with Mh ~ 125.7 and M^ ~ 98 GeV is certainly possible (see 
also [TUlIlS]), it is clearly not favoured by our rate analysis. 

The most relevant parameters for higher-order corrections in the MSSM Higgs sector are the soft 
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Figure 11: Stop mixing parameter X t /Mq 3 vs. the light stop mass (left), and the light vs. heavy stop 
masses (right) in the light Higgs case. 




Figure 12: Stop mixing parameter X t /Mq 3 vs. the light stop mass (left), and the light vs. heavy stop 
masses (right) in the heavy Higgs case. 

SUSY-breaking parameters in the stop sector. As it was shown in [10], light scalar top masses down to 
150 GeV are in agreement with a light CP-even Higgs mass around ^125 GeV, provided the mixing 
in the scalar top sector is sufficiently strong. Here we show the corresponding results including the 
rate analysis. In Fig. [ED we show X t /Mq 3 vs. the light stop mass (left plot) and the light vs. the heavy 
stop mass (right plot) in the light Higgs case. In the left plot one can see that the case of zero stop 
mixing in the MSSM is excluded by the observation of a light Higgs at ~ 125.7 GeV (unless Mq 3 
is extremely large), and that values of \X t /Mq 3 \ between ~ 1 and ~ 2.5 must be realised. For the 
most favoured region we find X t /Mq 3 = 2 — 2.5. It should be noted here that large values of \A t \ 
(\A t \ > ^/6Mq 3 ) could potentially lead to charge and colour breaking minima [57] . We checked that 
applying a cut at \A t \ > \^6Mq 3 would still leave most points of the favoured region. Concerning 
the value of the lightest scalar top mass, the overall smallest values are found at mi ~ 200 GeV, in 
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Figure 13: Dependence of the rates i?^ 7 and R bb (VH) on the stop mixing parameter X t /Mq 3 for the 
light Higgs case. 

agreement with [10] H Even taking the rate information into account, the (most) favoured values start 
a ^ ™>ii ~ ^00 GeV for positive X t . Such a light t\ is accompanied by a somewhat heavier £2, as can be 
seen in the right plot of Fig. [TTJ Still, values of ~ 200 GeV are realised for ~ 600 GeV, which 
would mean that both stop masses are rather light, offering interesting possibilities for the LHC. The 
highest favoured values we find are ^1.4 TeV. These are the maximal values reached in our scan, 
but from Fig. [TT] it appears plausible that the favoured region extends to larger values of both stop 
masses. Such a scenario would be extremely difficult to access at the LHC. For the intrepretation of 
these results it is important to remember that we have assumed a universal value for the soft mass 
parameters in the scalar top and bottom sector. Relaxing this assumption would potentially lead to 
larger regions of parameter space in which all applied constraints can be satisfied. 

The results for the scalar top masses in the heavy Higgs case look in principle similar to the 
light Higgs case, but with substantially smaller favoured regions, which are nearly solely realised for 
positive X t with X t /Mq 3 = 2-2.3, as can be seen in Fig. [T2l The favoured values of raj range between 
~ 250 GeV and ~ 700 GeV in this case, whereas the preferred range of the heavy stop extends from 
mi 2 rsj 650 GeV to ~ 1100 GeV. We now turn to the analysis of rates as a function of the 
underlying MSSM parameters. This comparison allows to analyse the various mechanisms that are 
responsible for the observed differences in the decay rates with respect to the SM values. In Fig. fT3l we 
analyse the correlation between the ratio X t /Mq 3 and R^ (left) or R bb (VH) (right) in the light Higgs 
case. It can be seen that the enhancement in the 77 channel is only substantial for Xf/Mq 3 > 2, where 
values of up to ~ 1.7 can be reached in the favoured region. Such an enhancement can have two 
sources: a suppression of T(h —> 66), as the by far largest contribution to the total width, or a direct 
enhancement of T(h —> 77). That the first mechanism is indeed responsible for a substantial part of 
the scenarios with an enhancement of R^ can be seen in the right plot of Fig. [T3l which together with 
the middle plot of Fig. [5] illustrates that the enhancement in the 77 channel in the favoured regions 
is accompanied by some suppression of the bb channel. This suppression/enhancement is realised for 
large values of X t /Mq 3 . 

5 In this work a Higgs mass measurement of Mh = 125 GeV was assumed, whereas we now use the average mass 
M H 125.7 GeV. 
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Figure 14: Correlation of the fi parameter to the value of Ma (left), and dependence of A b corrections 
on ^tan/3 (right), both in the light Higgs case. 



Such a suppression of the bb channel can happen in two different ways. In the MSSM, the (effective) 
coupling g hb - b is given by, 



9hbb 



1 + A 6 



sm% , a coso: eff 



cos /3 



sin/3 



(14) 



where a is the mixing angle in the CP-even Higgs sector. Higher-order contributions from Higgs 
propagator corrections can approximately be included via the introduction of an effective mixing 
angle, corresponding to the replacement a — >> a e ^ [39] (in our numerical analysis we treat propagator- 
type corrections of the external Higgs bosons in a more complete way, which is based on wave function 
normalisation factors that form the Z matrix [27] ). A suppression of the h — >> bb channel thus occurs 
for scenarios with small a e g. Furthermore, genuine corrections to the hbb vertex enter Eq. (TT4l) via 
the quantity A5 oc \i tan/3 [581 159] 1^1 see [11] for a more detailed discussion of the possible mechanisms 
giving rise to suppression of the h —> bb channel. 

While the loop-corrected coupling g hbb , as given in Eq. (Q3D, approaches the SM coupling in the 
decoupling limit (Ma ^> Mz), a suppression of g hbb is possible for Ma not too large provided that 
A b is numerically sizable and positive. We analyze this in Fig. [TU The left plot in this figure shows 
that the most favoured regions are obtained for fi > 1 TeV, and that the combination of small fi and 
Ma < 500 GeV is disfavoured. The corresponding A b values are shown in the right plot as a function 
of /xtan/3. The most favoured regions here are found in the range 0.3 < A b < 1.5, for correspondingly 
large values of /itan/3 ~ 30-70 TeV. Note that the large values for the A b corrections do not pose 
problems with perturbativity, since they tend to reduce the bottom Yukawa coupling. It should be 
noted that the A b corrections in Eq. (|T4l) have another important effect: while in the absence of 
those contributions a small value of would give rise to a simultaneous suppression of the Higgs 
couplings to bb and to r + r _ , the A b corrections differ from the corresponding contributions to the 
gh r + r - coupling. This implies in particular that the g^r+r- coupling may be suppressed while the g hbb 
coupling remains unsuppressed (and vice versa), see the discussion of Fig. [5] above. 



6 The dominant contributions to A5 beyond one-loop order are the QCD corrections, given in [60]. Those two- loop 
contributions are not included in our analysis, but their numerical effect is approximated by using a scale of m t for the 
evaluation of the one-loop expression. 
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Figure 15: Enhancement of the h — > 77 partial width in the presence of light staus for the light Higgs 
case (left) and heavy Higgs case (right). 

For the second mechanism, direct enhancement of the T(h — > 77) width, it is known that other 
SUSY particles can play an important role. One possibility that has been discussed recently is to have 
very light scalar taus [13]. The effect of light scalar taus can also be observed in our analysis, as can 
be seen in Fig. |T5j Here we show T(/i, H —> jj)/T(h, H —> 77)sm as a function of m^. In the light 
Higgs case, shown in the left plot, for ~ 100 GeV the enhancement over the SM width reaches 
50% in the favoured region. Even lower values of (which are allowed regarding the limits from 
direct searches, see [24]) are forbidden in our scan from the requirement that the LSP is the lightest 
neutralino, together with the lower limit of M2 > 200 GeV and the GUT relation between M\ and M2. 
Relaxing these assumptions would allow for a larger enhancement of T(h — > jj)/T(h 77)sm 5 as is 
clear from the sharp rise of this rate seen in Fig. [15] for low • For > 300 GeV a decoupling to the 
SM rate is observed. Through the contributions of light scalar taus it is thus possible to accommodate 
enhanced values of i?^ , while maintaining R% b and Ry V at the SM level (as also observed in our 
results, cf. Fig. [5|). In the heavy Higgs case, on the other hand, as shown in the right plot of Fig. [T5l 
the favoured region is located close to one, and light staus do not contribute to a possible enhancement 
of R* T 

Similarly to the light Higgs case, we investigate the dependence of the rates on the stop sector 
parameters for the heavy Higgs case. The results are shown in Fig. [161 As in Fig. [T2J, the favoured 
regions are given for large and positive X t /Mq 31 where we find 0.8 < RtL < 1.6 and a corresponding 
suppression of 0.6 < R^ b < 1.0. The A5 corrections, that enter analoguously to Eq. (fT4l) . can also in 
this case be largely responsible for the suppression of the R^ rate, as we show in Fig. [T71 Here one 
can see that in the heavy Higgs scenario only values of A5 between ~ 0.3 and ~ 0.6 are favoured, 
which are realised for 10 TeV < //tan/3 < 35 TeV, i.e. smaller values than in the light Higgs case. 

To summarize the discussion on favoured MSSM parameter regions, we list in Tab. [6]the parameter 
values for the best fit points (for the complete fit). We also give the parameter ranges corresponding 
to Ax^ H < 1. For several of the parameters this range extends to the limits of our scanned interval. 
Cases like this have been indicated in Tab. [6]with parentheses around the corresponding numbers. One 
can see that in most cases the ranges, even evaluated for Ax\ H < 1, are quite wide. One exception 
is tan/3, which is relatively tightly constrained (at least at the level of A%^ H < 1) in the light Higgs 
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Figure 16: Dependence of the rates RiL and Rff b (VH) on the stop mixing parameter X t /Mq 3 for the 
heavy Higgs case. 
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Figure 17: Correlation of the \i parameter to the value of Ma (left), and dependence of A5 corrections 
on //tan/3 (right), both in the heavy Higgs case. 



case, and even more so in the heavy Higgs case. In the latter case, as discussed above, also the masses 
of the additional Higgs bosons are relatively tightly constrained, offering some important information 
on how this scenario can be further explored (see the next section). More precise experimental data 
would be needed to achieve tighter constraints on the other fitted parameters, which enter the MSSM 
Higgs phenomenology via loop corrections. The fact that even in the more "exotic" scenario, where 
the signal at ~ 125.7 GeV is interpreted in terms of the heavier CP-even Higgs of the MSSM, the 
values of individual SUSY parameters are only moderately constrained by the fit illustrates that a 
reasonably good description of the data can be achieved without the need of tuning certain parameters 
to specific values. This is of course even more the case for the interpretation in terms of the light 
CP-even Higgs. 
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Lif 


rht Higgs 


case 


Heavy Higgs 


case 


Parameter 




Best fit 






Best fit 




M A [GeV] 


300 


669 


860 


120.5 


124.2 


128.0 


tan/3 


15 


16.5 


26 


9.7 


9.8 


10.8 


fi [GeV] 


1900 


2640 


(3000) 


1899 


2120 


2350 


Mq 3 [GeV] 


450 


1100 


(1500) 


580 


670 


740 


M~ h [GeV] 


250 


285 


(1500) 


(200) 


323 


(1500) 


Af [GeV] 


1100 


2569 


3600 


1450 


1668 


1840 


M 2 [GeV] 


(200) 


201 


450 


(200) 


304 


370 


M fe [GeV] 


122.2 


126.1 


127.1 


63.0 


65.3 


72.0 


M H [GeV] 


280 


665 


860 


123.9 


125.8 


126.4 


M^i [GeV] 


310 


673 


860 


136.5 


138.8 


141.5 



Table 6: Best fit parameter values (in the respective middle column) and ranges for H < 1. Values 
in parentheses indicate that the limit of the scan range has been reached. 

5 Conclusions and Outlook 

We have analyzed the compatibility of the Minimal Supersymmetric Standard Model (MSSM) with the 
recent discovery at the LHC of a Higgs-like state at Mh ~ 125.7 GeV. To this end we have studied the 
low-energy (phenomenological) pMSSM-7 parameter space, where we allowed the seven parameters 
most relevant for Higgs and flavour phenomenology to vary freely: the CP-odd Higgs boson mass, 
Ma, the ratio of the two vacuum expectation values, tan/3, a common soft SUSY-breaking parameter 
for the scalar top- and bottom quarks, M^ 3 , a soft SUSY-breaking parameter for the scalar tau and 
neutrino sector, M^, a common trilinear coupling for the third generation, Af, the higgsino mass 
parameter, /x, as well as the SU(2) gaugino mass parameter, M2. The U(l) gaugino mass parameter 
Mi was fixed from the value of M2 using the GUT relation. The other parameters have been set 
to fixed values as to be generically in agreement with recent SUSY searches at the LHC and with 
low-energy observables such as (g — 2) M . 

A random parameter scan over the seven free parameters with (9(10 7 ) scan points has been per- 
formed. For each scan point, a x 2 function was evaluated, taking into account the measured rates in 37 
individual Higgs search channels from ATLAS, CMS, and the Tevatron, the best-fit mass values of the 
LHC experiments, as well as the following low-energy observables: BR(£> — > X.7), BR(B S —> 
BR(B U rzv r ), (g - 2) M and M w . 

As a starting point we find that the SM yields a good fit to the data, with a y 2 per degree of 
freedom (dof) around unity. The precise value depends on whether low-energy observables and/or the 
Tevatron data are included in the fit. Turning to the MSSM, we find that the pMSSM-7 provides 
an excellent fit to the Higgs data in the case that the light CP-even Higgs is interpreted as the new 
state at ~ 125.7 GeV. In the case that the heavy CP-even Higgs boson is interpreted as the newly 
discovered state the fit is still acceptable, but somewhat worse than in the light Higgs case once 
Tevatron and low-energy data is included. The two MSSM best-fit points have a total x 2 /dof of 
30.4/36 (42.4/36) for the light (heavy) Higgs case, respectively, after the inclusion of LHC, Tevatron 
and low-energy data. This translates into p-values of 73% and 21%, respectively. The corresponding 
SM value for x 2 / l/ is 45.3/42, resulting in a fit probability of 34%. The largest x 2 contribution in 
the SM comes from the inclusion of (g — 2) M , which shows a more than 3 a deviation from the SM 
prediction. Regarding the comparison of the results for the light Higgs case and the heavy Higgs case 
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in the MSSM it should be noted that a sizable part of the additional x 2 contribution in the heavy 
Higgs case results from the BK(B U ru T ) measurement and from (g — 2) M . The agreement between 
theory and experiment (both for the MSSM and the SM) for BK(B U tv t ) would improve with 
the inclusion of the new Belle measurement. The x 2 contribution arising from (g — 2) M for the heavy 
Higgs case of the MSSM could potentially be improved if in addition to the seven parameters that are 
varied in our fit also the second generation slepton parameters would be treated as free fit parameters, 
which would essentially select the slepton mass parameters yielding the lowest x 2 value from (g — 2) M 
for each point in parameter space without affecting the other phenomenology. Thus, while the best 
description of the data is achieved if the new state at ~ 125.7 GeV is interpreted as the light CP-even 
Higgs boson of the MSSM, the more "exotic" interpretation in terms of the heavier CP-even Higgs 
of the MSSM is also permitted by the data, even if the results from the Higgs searches at the LHC 
are supplemented with results from the Tevatron Higgs searches and with results from flavour physics 
and electroweak precision data. The latter interpretation would imply that also the other four Higgs 
bosons of the MSSM would be rather light, giving rise to exciting prospects for the searches for non 
SM-like Higgses. 

In the case of the light CP-even Higgs with ~ 125.7 GeV, we find for the best-fit point in the 
full fit an enhancement of production times branching ratio for the 77 channels of about 40% with 
respect to the SM prediction. The rates for the gauge boson channels that we obtain are similar as 
in the SM, and the same holds for the fermionic channels (bb and t + t~). While the fit results for the 
77, VV , and bb rates show a clear x 2 minimum, the t + t - channel has a very broad distribution close 
to the minimum, and no strong preference can be attributed to the actual best-fit value. 

In the case of the heavy CP-even Higgs with Mh ~ 125.7 GeV we find for the best fit point a 
somewhat smaller enhancement of the 77 channel, an enhancement of the gauge boson channels, and 
a suppression of the t + t~ (VBF) and bb channels, whereas the t + t~ inclusive channel is enhanced 
due to the contribution of the CP-odd Higgs boson. 

For the light Higgs case, as well as for the heavy Higgs case, the rates in the 77 and VV channels are 
strongly correlated, however in most cases with the possibility of a stronger enhancement (or smaller 
suppression) in the 77 channel. Between the 77 channel and the bb channel an anti-correlation can 
be observed. This shows that the fit within the MSSM favours at least over a part of the preferred 
region a scenario where a 77 enhancement is caused by a suppression of the bb channel. Our results 
furthermore show that the t + t - channel can be strongly suppressed, while the bb channel can remain 
close to the SM strength. Since we find that the gg — >► h,H production channels in the favoured 
regions are not substantially enhanced, the observed enhancement of the 77 channel is not caused by 
larger production cross sections. A suppression of the bb channel in the light and the heavy Higgs case 
can be caused by a large value of A5, which can reach values exceeding unity in the favoured region. 
In the light Higgs case, the 77 channel can also be enhanced by the contribution of light scalar taus 
to the decay process. In the case where the lightest scalar tau mass is as low as about 100 GeV, we 
find an enhancement of up to 50% from this mechanism. 

For the scalar top masses, we find that the favoured regions start at ~ 200 GeV and ~ 
600 GeV in the light Higgs case. They extend up to ~ 1.4 TeV and ~ 1.6 TeV, respectively, which 
are the maximal values accessible in our scan. The mixing in the scalar top sector must exceed 
\X t /Mq 3 \ ~ 1, where the most favoured regions have X t /Mq 3 = 2 — 2.5. Similar values for the lower 
bounds on the scalar top masses are found in the heavy Higgs case. However, for this case we find 
that the favoured regions are also bounded from above by (roughly) m^ 2 < 1 TeV. 

As is evident from our analysis (as demonstrated e.g. by Figs. [1] and [2j), the fitted rates in the MSSM 
interpretations are not significantly different from the SM predictions, using the current experimental 
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and theoretical uncertainties. Therefore, if no other new states beyond the current Higgs candidate are 
discovered in the near future, the question arises how much the precision of the current measurements 
would need to be improved in order to distinguish the MSSM from the SM, based on precisely measured 
rates alone. In order to obtain a rough estimate for the answer to this question, we set the hypothetical 
future central values of the measurements to the MSSM best fit point in the light Higgs case, so that we 
can investigate the impact of prospective future experimental results in a scenario where this particular 
MSSM point is actually realized in nature. Then, we scale the full uncertainties of the signal rates 
measured in the LHC8 (i.e., the channels measured at the LHC with = 8 TeV) channels by a 
global scale factor and infer the value of the scale factor at which the deviation of the SM prediction 
from the assumed future measurements with higher precision reaches a significant level. We scale 
only the LHC8 results (and not LHC7, Tevatron, or LEO) to account for those channels for which an 
actual improvement of precision is expected in the future (at = 8 TeV or higher). To reach 2 a 
(3 a) significance for rejecting the SM, we find that a scale factor of 21 % (18%) is required. Reducing 
the uncertainties to ~ 20% of their current values for all channels (i.e., a factor five improvement 
compared to the present situation) appears to be reasonable as a long-term goal for the LHC [61J. For 
the necessary improvement of the statistical, systematical and theoretical precision to go significantly 
beyond this number, as needed in this case to go towards 5 a sensitivity for SM exclusion, would 
require something close to the expected experimental sensitivity of the International Linear Collider 
(ILC), see [62] and references therein. The corresponding numbers in the heavy Higgs case are 31% 
(26%), which would be needed to reach a 2 a (3 a) rejection of the SM via Higgs rate measurements 
alone. 

While distinguishing the fitted MSSM interpretations from the SM by virtue of improving the 
precision of the rate measurements alone may be difficult, promising paths to establishing the presence 
of physics beyond the SM in the Higgs sector are on the one hand the search for additional (non 
SM-like) Higgs states and on the other hand the precise measurement of further Higgs properties. 
Concerning the latter, a promising example would be to measure the Higgs CV properties in the 
t + t - final state (assuming that the existence of this decay mode will be confirmed for the observed 
signal). In the heavy Higgs case, as mentioned in Sect. 13. 2\ the contributions from H and A add up 
in channels with poor mass resolution when no qhvv coupling is involved. This is the case for the 
inclusive production mode with decays into t + t~. There, an analysis of the r spin correlations [63] 
(and references therein) could show significant deviations from the SM prediction of CV = +1, due to 
the presence of the A production mode. In contrast, in the light Higgs case, the discovered state at 
Mh ~ 125.7 GeV would be measured with exact CV = +1. 

The potential discovery of additional Higgs like states would be the clearest way to distinguish 
the MSSM h and H interpretations from each other. As discussed above, see in particular Tab. [6J 
the masses of the yet undiscovered Higgs bosons are expected in different ranges for the h and H 
interpretations. However, discovering the lighter CP-even state h in the heavy Higgs case appears to 
be difficult at the LHC due to its low mass and heavily suppressed coupling to vector bosons, g\yy. 
At the ILC, however, the hA and H + H~ production modes would be straightforward to measure in 
the heavy Higgs case. On the other hand, in the light Higgs case all other states could be beyond 
the mass reach within the current fit uncertainties. Here, the ILC measurements would have to rely 
on precise coupling measurements to distinguish the h and H interpretations directly (see [62J and 
references therein). In that case also a combined interpretation of LHC and ILC data might be very 
valuable [SUES]. 

In the heavy Higgs interpretation, the LHC searches for MSSM Higgs bosons in the r + r~ final 
states will soon have sensitivity to start probing the region of small Ma and moderate tan /3 that is 



28 



favoured in this scenario. Furthermore we find that the charged Higgs boson in this case should be 
lighter than mt, so that at the production from top quark decays at the LHC would be kinematically 
possible. The favoured values of BR(t — >> H + b) are just below the current experimental bounds. Still, 
an improvement of the limits by one order of magnitude would be required to fully cover this possibility. 

New data from the ATLAS and CMS Higgs boson searches is rapidly emerging. It will be particu- 
larly important to investigate on the one hand potential deviations of the rates from the SM predictions 
and on the other hand the outcome of searches for additional non SM-like Higgses. Confronting these 
results with predictions in the MSSM will show whether this model, whose unambigous prediction of 
a light (and potentially SM-like) Higgs boson seems to be well supported by the data, will continue 
to provide a viable description of nature also in the future. 
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